Granular discorect angle in a thermalized bath of hard disks 
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By using tiie Enskog-Boltzmann approacii, we study tlie steady-state dynamics of a granular 
discorect angle placed in a two-dimensional bath of thermalized hard disks. Hard core collisions 
are assumed elastic between disks and inelastic between the discorect angle and the disks, with a 
normal restitution coefficient a < 1. Assuming a Gaussian ansatz for the probability distribution 
functions, we obtain analytical expressions for the granular temperatures. We show the absence of 
equipartition and investigate both the role of the anisotropy of the discorect angle and of the relative 
ratio of the bath particles to the linear sizes of the discorectangle. In addition, we investigate a 
model of a discorectangle with two normal restitution coefficients for collisions along the straight 
and curved surfaces of the discorectangle. In this case one observes equipartition for a non trivial 
ratio of normal restitution coefficients. 

PACS numbers: 05.20.-y,51.10+y,44.90+c 
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I. INTRODUCTION 

Granular matter is characterized by the existence of 
dissipative forces between particles. In order to sustain a 
collective motion, it is necessary to provide energy con- 
tinuously. When power supply is sufficiently copious, the 
assembly of granular particles attains a non equilibrium 
steady state (NESS)[l|, I3- It is customary to define the 
granular temperature as the second moment of the ve- 
locity distribution. It is a source of both fascination and 
inconvenience that the well-know properties of temper- 
atures characterizing thermal systems are not necessar- 
ily transferable to granular temperatures. In particu- 
lar, recent work both theoreticalU HHH S lE 13 and 
experimental |9l Uojl . has shown that in a binary granular 
system the two species have different granular temper- 
atures that are non-trivial functions of the microscopic 
parameters (mass, size, restitution coefficient,...). Al- 
though the absence of equipartition is not surprising for a 
dissipative system sustained in a NESS, a more complete 
investigation is necessary since the granular temperatures 
play an important role in hydrodynamic descriptions of 
these systems (in particular the absence of equipartition 
in binary mixtures yields granular temperature gradients 
which enhance segregation llj). In addition, the exten- 
sion of the fluctutation-dissipation theorem is an impor- 
tant issue in the context of granular gases^^- Other 
consequences of the absence of equipartition include the 
ability of a binary system to exhibit a segregation phe- 
nomena in a "Maxwell demon" experiment 13]. See also 
the homogeneous cooling state of a granular mixture |14| 
and the impurity problem[l5j. 

Most of the above-referenced studies examined assem- 
blies of spherical particles. Yet, in reality, the par- 
ticles composing granular systems are to some degree 
anisotropic and, in many cases, strongly so. Even if the 
particles are smooth, each collision results in some ex- 
change and, possibly loss, of rotational kinetic energy. 
There are relatively few studies of these systems (but see 



for example, Refs.[Tl [H El) and fewer still that focus 
specifically on equipartition. Huthmann et al. used 
kinetic theory to examine the free cooling of a system of 
granular needles in three dimensions and, more recently, 
two of the present authors studied a two-dimensional sys- 
tem composed of a single granular needle in a thermal- 
ized bath of point particles 20] in a NESS. For inelastic 
needle-point collisions, the rotational granular tempera- 
ture is smaller than the translational one while both are 
less than the bath temperature. The validity of the the- 
oretical predictions were confirmed by comparison with 
numerical simulations of the model. While this study 
provided useful insights, infinitesimal width of the parti- 
cle is obviously an idealization. 

The objective of present article is to consider a more 
realistic system where both the tracer particle and bath 
particles are of finite extent. Specifically, we consider a 
discorectangle in a bath of thermalized hard disks. For- 
tunately, despite the increased complexity, it is still pos- 
sible to obtain an analytic solution of the steady state 
kinetic equations. The principal difference between the 
discorectangle-disk and needle-point systems is that two 
kinds of collision are possible in the former compared to 
one in the latter. Specifically, a disk can collide with ei- 
ther the sides or the caps of the discorectangle. If each 
type of collision is characterized by different normal resti- 
tution coefficients we show that equipartition between 
the translational and rotational degrees of freedom can 
be obtained for specific values of these parameters. Con- 
sequently, for appropriate ranges of the restitution coeffi- 
cients the translational granular temperature may be less 
than or greater than the rotational one. 



II. MODEL AND COLLISION RULES 

We investigate a two-dimensional system consisting of 
a discorectangle of total length L -\- 2R, radius R and 
mass M with a moment of inertia / (The value of which 
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is given in appendix A). The bath consists of disks of 
mass m and of radius r. The vector positions of the 
center of mass of the discorectangle and a disk particle 
are denoted by ri and r2, respectively. The orientation 
of the discorectangle is specified by a unit vector Ui that 
points along the long axis. Let ri2 = ri — r2 and 
denote a vector perpendicular to ui . A collision between 
a discorectangle and a disk can take place either on the 
linear part or on the circular parts of the former, 



(a) 



ri2.u^ = -{R 



if |A| < L/2, and 



ri2-Ur 



-{R + r), 



(1) 



(2) 



if R > |A| — L/2 > 0, where denotes the unit vector 
of the collision axis, (see Fig.^. The relative velocity of 
the point of contact V is given by 



V = Vi2 + LJi X OC, 



(3) 



where Wi denotes the angular time derivative and OC 
the vector from the center of the discorectangle to the 
point of impact. 

The pre- and post-coUisional quantities (the latter are 
labeled with a prime) satisfy: 



Total momentum conservation 



Afv'l + r7iv'2 — A/Vl + TOV2. 



(4) 



• Angular momentum conservation with respect to 
the point of contact 

Ioj[k = lujik + mOC X (v2 - v'2) , (5) 

where k is a unit vector perpendicular to the plane. 

As a result of the collision, the relative velocity of the 
contacting points changes instantaneously according to 
the following relations: 



V'.u^ 
V'.ui 



(6) 
(7) 



where a is the normal restitution coefficient and u;j'- and 
Ui denote the unit vectors of the collision along the rec- 
tilinear part of the discorectangle. When the collision 
occurs on the circular parts of the discorectangle, the 
collision rules are given by 



V'.Ur = —aV.Ur, 
V'.ue = V.ue, 



(8) 
(9) 



where and ug denote the unit vectors associated with 
the circular part of the discorectangle (See Fig^jD). 

The tangential restitution coefhcient is set to one for 
the sake of simplicity. This choice is reflected in the form 
of Eqs. (EEJ. 
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FIG. 1: Geometry of the discorectangle and a disk in the 
plane: (a) Collision between the rectilinear part of the dis- 
corectangle and a disk, ri2 denotes a vector joining the point 
labeled 2 and the center of the discorectangle; ui is a unit vec- 
tor along the axis of the discorectangle, A is the algebraic dis- 
tance between the center of the discorectangle and the point 
of impact and Ui is a unit vector perpendicular to the axis of 
the discorectangle. (b) Collision between the circular part of 
the discorectangle and a disk. 9 denotes the angle between the 
direction of the long axis of the discorectangle and the axis 
defined by the contact point and the center of the disk.Ur and 
ug are unit vectors normal and tangential to the surface at 
the point of contact, respectively. 



By combining Eqs. lj2Jl- ® one obtains, after some 
algebra, the change of the discorectangle momentum 
Ap — M{v[ — Vi) for a collision along the linear part 



^_ (l + a)V.UjL 



rn ' M ^ / 



(10) 



for |A| < L/2 and at the two ends of the discorectangle 

(11) 



(1 + a)V.Ur 



M ' 4/ 



for i? > |A| - L/2 > with cos(6l) = 
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III. BOLTZMANN EQUATION 

Since we are interested in the homogeneous state, 
the distribution function /(vi,u;i) of the discorectangle 
obeys 



df{vi,uji) 
dt 



N / ^ / d-v-2 / dr2ri2/(vi, wi, V2) 



(12) 

where N is the total number of disks, /(vi, wi, V2) is the 
distribution function of the discorectangle and a disk, and 
T12 is the collision operator between a discorectangle and 
a disk. 

Defining the granular temperatures as quadratic av- 
erage of the appropriate velocity distribution, one has 
Tt = M/2(v2) and Tr = /(cj^) for the translational and 
rotational granular temperatures, respectively (the angu- 
lar brackets denote the average). By taking the second 
moment with respect of the velocity and of the angular 
velocity of Eq. H12() , one obtains 

2dTT 



Mdt 

= N 



dvidtJiat(vi/(vi, wi)) 



2 / dr2ri2/(vi,u;i, V2)v^, 
(13) 



dt 

|A|)e(|A| 

d\ri2.Ur\ 



dt 



where 612 is an operator that converts pre-coUisional to 
post-collisional quantities. 

The others terms of the collision operator corre- 
spond to the necessary conditions of contact Q{L/2 — 

|A|)(5(|ri2.Uj^| — (r + R)), and approach 9 ^ 

in the first case and Q{R + L/2 — 

L/2)(5(|ri2.Ur| — {r + R)), and approach Q 

in the second case. 

By taking the second moments of the distribution func- 
tion of the discorectangle and after substitution of the 
collision operator (Eq.l^SJ), one obtains explicitly for the 
translational kinetic energy 



J ...J dr2dvidv2dcjid0i [e(L/2 - |A|),5(|ri2.u^| - (r + R)) 



d|ri2.UjL 



dt 



e - 



d|ri2.uf 



dt 



/(vi,cji)$(v2)A^; 



+ e(i? + L/2 - |A|)e(|A| - L/2),5(|ri2.u,| - (r + R)) 



rf|ri2.Ur 



dt 
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d|ri2.Ur 



dt 



/(vi,cji)$(v2)A£;; 



0. 

(17) 



= j d-vidwidt{ujlf{vi,uJi)) 



Idt 

N I dvi I duJi 



d9i 
2tt 



dv2 I dr2Ti2/(vi,ti;i, V2)tJi. 

(14) 



In the stationary state the time derivatives of the left- 
hand side of these two equations are equal to zero. 

To build the collision operator between the discorect- 
angle and a disk, T12, one must include the change in 
quantities (i.e. velocity and angular momentum) pro- 
duced during the infinitesimal time interval of the colli- 
sion. This operator is different from zero only if the two 
particles arc in contact and if the particles were approach- 
ing just before the collision 21j. For a collision between 
a disk and the rectilinear part of the discorectangle, the 
explicit form of the operator given by 



ri2 «e(L/2- |A|),5(|ri2.u^| -r 



d|ri2.u^ 



dt 



e - 



d\ri2.uj^ 



dt 



R) 
(612 



1), (15) 



where &12 is an operator that changes post-collisional 
quantities to pre-coUisional quantities and Q{x) is the 
Heaviside function, and for a collision at two ends of the 
discorectangle 

T12 cx e(i? + L/2 ~ |A|)e(|A| - L/2),S(|ri2.u,| - (r + R)) 



d|ri2.Ur 



dt 



Q - 



d|ri2-Ur 



dt 



(&12-I), 



(16) 



A similar equation can be written for rotational kinetic 
energy. Since the impulse of the collision depends on the 
location of the impact, it is easy to show that the solu- 
tion is not a Maxwell distribution function, a property 
already observed in the model of a needle and points|20j. 
However, since the deviations from the Maxwell distribu- 
tion are small, we use it as a trial function. That is we 
assume that 



/(vi,u;i) cx exp 



IluI 



2^tT 2^rT 



(18) 



where 7^ and are the ratios of the translational and 
rotational discorectangle temperatures to the bath tem- 
perature, respectively. In summary, in order to obtain 
the granular temperatures of the discorectangle, it is nec- 
essary to: (i) Calculate the change of translational and 
rotational energy occurring during a collision (ii) Perform 
an average over all degrees of freedom of the collision in- 
tegral. 



IV. CALCULATION AND RESULTS 

A. Energy changes during a collision 

When a disk collides with the discorectangle, the 
change of the translational kinetic energy of the latter 
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is given by 



Ap.vi 



1 Ap2 



M 2 



-il + a)- 



for lAI < and 



V.ufvi.u^L 1 (l + a)2(V.uf 



±\2 



J_ 

M 



X2_ 
I 



2M 



j_ 

M 



A2N 



(19) 



(1 + a)V.UrVi.Ur 



A/ 



41 



+ 



1 (l + a)nV.u 

V m ^ M ^ 4/ / 



(20) 



for i? > |A| - f > 0. 

The collision results in a change of rotational energy, 
for lAI < ^, 




A(l + a) V.UjL(cj; +wi) 



M 



hi 
I 



-X{l + a)- 



V.UjLt^l A2(l + a)2(V.UjL)2 



^ + J_ + ^ 

™ ^ M ^ / 



Vm ^ M ^ / 



2\2 



(21) 



and for i? > |A| - L/2 > 0, 
ASf = -(1 



Lsin 6'V.u, 



2 

\2r2„;„2, 



A/ 



L2f 



41 



{l+afL^ sin^6l(V.u^ 



8/ 



\ m 



J_ 

M 



41 



(22) 



B. Expressions of the granular temperatures 

After inserting Eq. p8|) in Eq. (|17|l it is necessary 
to perform integrations over each variable in the corre- 
sponding equation for the rotational energy. Details of 
this rather technical task are given in the Appendix B. 
After some tedious calculation, one obtains the following 
set of equations 
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FIG. 2: Ratio of the translational (full curve) 7t and rota- 
tional 7_R (dashed curve) granular temperature to the temper- 
ature of the bath versus the normal restitution coefficient q 
for a system of a needle and a bath of points and a system of 
a needle and a bath of disks for c = 2/3 (a) and for c = 1/10. 
The curves corresponding to the needle-point system are al- 
ways above those for the needle-disk system (b). The insets 
display the differences 7t — 7i? versus a for the needle in the 
bath of disks (full curve) and the needle in the bath of points 
(dashed curve). 



6(c/i°i(a,/fc)-f (l-c)jf (a,/fc)) = 

i±^(c/f(a,fc) + (l-c)jf(a,fc)) 
a (c/i" (a, /c) + (1 - c) Ji" (a, k)) = 



(23) 
(24) 



where 



fc = 



M + m 



M -f m7T ' 



(25) 
(26) 
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M + m 



It 



M - 



and 



L/2 
L/2 +ir + R) 



and and are given by 



I'^iu.v) = dx- ' 
Jo 



(1 



v,2\rn 



Jo (1 + t^sm 0)™ 



(27) 
(28) 

(29) 
(30) 



Explicit expressions for the integrals appearing in 
Eqs. (Eni) and are given in Appendix C. Eq. 
is an implicit equation for a that, for a given value of 
a, can be solved with standard numerical methods, b is 
then easily obtained by calculating the ratio of integrals 
of Eq. (|23|) . Finally, from the values of a and 6, 7t and 
7_R, can be obtained from Eas. l|^ - l(77|l . 

A first check is the elastic case, a = 1, for which 
one obtains a = 6 = 1 which gives 7^ = 1 and, since 
a/h = ^r/^t, 1r — ^ i-e., the temperatures of trans- 
lational and rotational degrees of freedom are the same 
and correspond to the bath temperature, a property of 
an equilibrium system. 

In the limit i? — * and r — > 0, for which c ^ 1 corre- 
sponding to a needle in a bath of point particles, one re- 
covers the results of Rcf 20] . More interesting is the limit 
i? — > with r remaining finite, i.e., a needle in a bath of 
disks. By using Ea. (|28|l one obtains that c = L/{L + 2r). 
Since c < 1, unlike the simple needle-point system, there 
is a contribution resulting from collisions between bath 
particles and the needle's tips. This type of collision is 
dominant when the bath particles are larger than the 
longest size of the anisotropic particle. 

To illustrate the effect of this contribution. Figures |2K 
and b compare the system of a needle in a bath of point 
particles and the system of a needle in a bath of disks, the 
radius of disks r are equal to L/4 in Fig|2jL and r = 9L/2 
in Fig|23. It is noticeable that the translational and rota- 
tional temperatures of the latter system are smaller than 
those of the needle in a bath of points, an effect that 
becomes more pronounced when the radius of the disks 
becomes larger than the length of the needle. The dif- 
ference 7t — 7_R are shown in the insets for the needle 
and the bath of points and for the needle and the bath 
of disks. Note that the translational temperature is al- 
ways larger than the rotational temperature and that the 
difference depends very weakly on the size of the radius 
of the bath. Additionally, the difference increases when 
the restitution coefficient decreases from 1 (elastic case), 
reaches a maximum for a value of a ~ 0.3 and decreases 
slightly when the restitution coefficient still decreases. 

Figure 13 shows the influence of the anisotropy of the 
tracer particle. The rotational temperature is still lower 
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FIG. 3: Ratio of the translational (full curve) 7t and rota- 
tional 7ij (dashed curve) granular temperature to the tem- 
perature of the bath versus the normal restitution coefficient 
a for a discorect angle and a bath of disks, for different values 
of the anisotropy r = R and L = R, 3R, lOR 



than the translational temperature what ever the elon- 
gation of the particles, but the effect is greater when the 
elongation is large. The two upper curves correspond the 
translational (full curve) and rotational (dashed curve) 
temperature of the discorectangle when L = R, the 
two intermediate curve to granular temperatures when 
L ~ 3R, and the two lower curves for a discorectangle 
with L = lOi?. 

When the anisotropy of the discorectangle approaches 
zero, i.e. L/R —fO,c—>0, one can show by using 
Eqs.^-^ that 



I 



7t 



2+t(l-")' 



(31) 



which is the the result of Martin and FiaseckiQ for a 
spherical tracer particle in a bath of spherical particles. 
Moreover, the limit L/R — > leads equipartition between 
the rotational and translational granular temperatures, a 
result which is different from a model of pure spherical 
particles where, in the absence of tangential friction, the 
particles cannot exchange rotational energy. In physi- 
cal systems, the particles are never completely spheri- 
cal and our model shows that if an infinitesimal amount 
of anisotropy is present, the translational and rotational 
temperatures are equal in the steady state. Although our 
analysis provides no quantitative information about the 
relaxation time to reach this NESS, it is certain to be 
very long in this limit of small anisotropy. 
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FIG. 4: Ratio of the granular temperatures 7t (full curve) 
7_R (dashed curve) to the bath temperature as a function of 
the mass ratio m/M , for a homogeneous discorectangle with 
c = 4/7. 

C. Influence of mass ratio 

We consider a homogeneous discorectangle of mass M 
in a bath of disks each of mass ni for which AI ^ m. 
When m/M one obtains, from Eqs.|2Sl)-(|21|l, that 

1 + a 

1T = 1R^ (32) 

i.e., equipartition between the degrees of freedom of the 
tracer particle, but not between the bath and the tracer 
particle. This is, moreover, the same result for a needle 
in a bath of point particles. We conjecture that this 
result is general in the sense that we expect equipartition 
between the different degrees of freedom of the tracer 
particle in a bath of light particles, whatever the shape 
of the tracer particle and the dimension of the system. 
The behavior for finite values of the ratio m/M is shown 
in Fig^ The granular temperatures decrease when the 
ratio m/M increases, and the translational temperature 
remains higher than the rotational temperature for each 
value of a. 



V. NON UNIFORM RESTITUTION 
COEFFICIENT 

In practice it may be difficult to construct a discorect- 
angle for which the restitution coefficient is constant over 
the entire perimeter. It is clear, for example, that if the 



FIG. 5: Sketch of the discorectangle with two restitution co- 
efficients ai (rectilinear part) and ct2 (circular part). 
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FIG. 6: Ratio of the translational (full curve) 7t and rota- 
tional 7_R (dashed curve) granular temperature to the temper- 
ature of the bath versus the normal restitution coefficient ol2 
with a normal restitution coefficient a.\ — 0.5 for a homoge- 
neous discorectangle with M — m, L — 2R and R = r. Note 
that equipartition is recovered for a non-trivial value of 02 



object is composed of a homogeneous viscoelastic ma- 
terial, collisions with the ends will be characterized by 
a smaller restitution coefficient than collisions with the 
linear part. This effect becomes more pronounced as the 
elongation {L/R) increases. Other possibilities exist for 
a non-homogeneous discorectangle composed of two or 
more materials. For example, a hard material may be 
used to construct the caps. In addition, the restitution 
coefficient could depend on the relative velocity of the 
point of impact ^231, an effect that one neglects here as a 
first approximation. 

As a first approach to describe this possibility, we con- 
sider in this section a discorectangle where the restitution 
coefficient is equal to ai for a collision along the rectilin- 
ear part of the object and equal to a2 for one along the 
circular part: see Fig[Sl Using the procedure outlined 
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above one obtains the following set of closed equations 
b ((1 + ai)cl^\a, k) + {l + a2)(l - c)J^\a, k)) 

(33) 

a ((1 + ai)cll\a, k) + {1 + a2){l - c)Jl\a, k)) 

- ^^^cll^a^k) + - c)JlHa.k) 

(34) 

Figure El shows the ratio of translational and rotational 
temperature of a discorectangle with R — r and L — 2R 
to the temperature of the bath as a function of the nor- 
mal restitution coefficient a2 with a fixed ai = 0.5. One 
notes that the translational temperature becomes smaller 
than the rotational temperature for a2 > 0.765. When 
the two curves cross, equipartition is recovered, but un- 
like the limiting cases discussed above of a discorectangle 
with a uniform restitution coefficient (light bath parti- 
cles, infinitely small anisotropy), for a non-trivial value 
of restitution coefficient a2 and, moreover, for larger val- 
ues of a2 the ratio of temperatures is inverted. 

One can determine in general when equipartition is 
recovered for a discorectangle with two restitution coeffi- 
cients. Using the relation a = b (assumption of equipar- 
tition) and Eqs.(j23)- one obtains two implicit equa- 
tions with the three parameters ai, a2 and a. A simple 
numerical procedure allows us to obtain the a2 as a func- 
tion of ai. 

Figure 13 shows the equipartition lines in the (ai,a2) 
space. Above each line, corresponding to a given elonga- 
tion of the rectangle, Tr > Tt while the reverse inequal- 
ity applies to the region below the line, ft is noticeable 
that as the elongation increases, the region of the (ai , 02) 
space where Tji > Tt decreases. 

Figure IHI shows the role of the size of the bath particles 
on the existence of equipartition of a discorectangle of 
length L = 8R. For small bath particles, only a small 
range of ai, ( between and ~ 0.3) with a smaller range 
of a2 ( between ~ 0.89 and 1) allows equipartition for the 
tracer particle. For larger bath disks, all values of ai (be- 
tween and 1) are available with a smaller corresponding 
range of a2. 

VI. CONCLUSION 

We have investigated the influence of the anisotropy 
of a tracer particle in a bath of thermalized disks in two 
dimensions. By using a mean- field approach, we have 
obtained analytical results for the rotational and trans- 
lational temperatures. For a homogeneous discorectangle 
with a uniform normal restitution coefficient, the trans- 
lational temperature is always higher than the rotational 
temperature, with the difference depending on the elon- 
gation of the tracer particle, the size of the bath particles 
and the mass ratio. 




FIG. 7: ai versus 02 where equipartition is obtained for r = R 
and different values of L = R, 2R, . . . , 5R from bottom to top 



At present, the exact dependence of the restitution co- 
efficient on the position of the point of impact is un- 
known, ft seems likely, however, that it is not constant 
along the perimeter. As a first approximation, we con- 
sidered an extension of the model where it takes on two 
values, one for collisions with the linear part and another 
for collisions with the caps. In this case the difference 
between the translational and rotational granular tem- 
peratures can be either positive or negative depending 
on the values of the restitution coefficients and the other 
parameters. It is clearly possible to generalize the cal- 
culation to allow for a continuously varying restitution 
coefficient. 

It would also be interesting to investigate a system 
of a tracer particle with a small anisotropy where the 
tangential restitution coefficient has a non-trivial value 
(— 1 < aT < 1) where the limits correspond to a perfectly 
rough and a perfectly smooth surface, respectively. The 
intermediate situation corresponds to some friction which 
is important in some granular systems. 



APPENDIX A: MOMENT OF INERTIA 

For a homogeneous discorectangle, the moment of in- 
ertia is given by 

lo.^ J J^p{x^ +y^)dxdy (Al) 
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FIG. 8: a\ versus 02 where equipartition is obtained for L — 
8R and for different values of disk radius r/R = 0,1,2,3,5 
(number labelling each curve) 



and the mass of the system is 



M 

which gives 



pdxdy = p{ttR^ + 2LR). (A2) 



loz = pR 



ttR 



i\2 



-m^ (A3) 



By substituting the density as a function of the total 
mass of the discorectangle. 



I ^ M 



3 



vri?(f + 



TrR + 2L 



(A4) 



Two well-known limits are recovered 



lim /, 



Oz 



lim loz = 



12 



(A5) 
(A6) 



For an inhomogeneous discorectangle, the moment of in- 
ertia depends on the mass distribution. However, it is 
possible to determine the lower and upper bounds for al- 
lowable values. A trivial lower bound of zero is obtained 
when the mass is concentrated at the center. Conversely, 
when the mass is distributed equally at the two extrem- 
ities of the object (point masses of Af/2 at a distance of 



L + R from the center on each side) , one obtains 

I = M{LI2 + Rf (A7) 
which gives the upper bound for the moment of inertia. 



APPENDIX B: NEEDLE AVERAGE ENERGY 
LOSS 

As for binary mixtures of spheres Q, we use a Gaussian 
ansatz for the distribution functions and introduce two 
different temperatures corresponding to the translational 
and rotational degrees of freedom of the needle. The 
homogeneous distribution functions of the needle and of 
the points are then given respectively by 



/(vi,cji) exp 



2T 



2T 



$(v2) ~ exp 



2T 



(Bl) 



(B2) 



where T is the temperature of the bath, 77^ and 7/? the 
ratio of the translational (and rotational) temperature of 
the needle to the bath temperature. 

We introduce the vectors x and v such that 



1 



^2T{M-fT + m) 
I mM 



(Mvi + TOV2) 



(vi - 7TV2) 



(B3) 



(B4) 



The scalar products V.uj'- and V.u^ can be expressed as 



V.UjL = h 
V.u,. — h 



(7T - i)x-ui + Vrr 

A 

(7T " l)x.Ur + + \j ^)'"-"' 



-^1 

(B5) 



' oji — sm f 



where h = 



2T 



A/7T +m ' 



Let us introduce ^ = oJi 



(B6) 



The translational 



energy loss is given by the formula 



p=±i 



\V.ut\e{pV.nt)ei- - \X\)AEt 



{R + r) / d0 



dOi 
"2^ 



V.Ur\e{pV.Ur)AE{ 



dx i dv I d^ exp(-x ~ ly - ^ 



(B7) 



= 



9 



Since Eq. ljB5|) depends only on X-'^i' a-nd v.uj- , one 
can freely integrate over the direction of Ui for the vectors 
X and v, and similarly for Ea. (|B6(l . The integration over 
9i can be easily performed. If we introduce the three 
dimensional vectors G^j- and s^j_ with components: 



APPENDIX C: INTEGRALS 



The coupled equations (|^ - 
eight integrals /?"^(a, k), /^'^(a, k), 



depend on 
{a,k),I^^a,k) 



the 



2T 



{it - 1), 



and 



(B8) 



J°^(a, fc), J2^(a, k), Jl^{a, k), J2^{a, k) which can be ex- 
pressed in terms of transcendental and special functions. 
For completeness, we give below their expressions. 



M-fT + m y/^W' V / J Ii\a,k) = 
(B9) 



ln{Vak + \/l + ak) 



1-a / /(l-a)fc 

— - — arctan \ — — 

k \ V 1 + afc 



(CI) 



= (X.UjL,i/.UjL,0 



(BIO) 
(Bll) 



Respectively, one has introduces Gu^ and vectors 
associated with collisions on the circular parts of the dis- 
corectangle by changing in and A by ^sin0. By 



inserting Eq. lO in Eq. ((B7)| . 

can be rewritten as 



the average energy loss 



,3/2 



^/k 



2VW 



: arctan 



I^^{a, k) = —j=- ln(Vafc + VT+ak) + 



(1 - a)^/TTak 



2{l + k) 



l{l-a)k 
1 + afc 



(C2) 



1 (l + a)2(G„^.s, ,)2 



/i (a,fc) 



2M (J--)- 1+^)2 
(l + a)G„i.s„i . 



2T 



-ilTSl 



M 



-S2 



+ {R + r) de J ds„^cxp(-s2j|G |e(pG„^.s„^) 
1 (l + a)2(G„„.s„„)2 



2M fJ- + J_ + iisini8\2 
V m ^ M ' 4/ 

(l + a)G„^.s„^ / 



2T 



m ^ Af ^ 4? V ' 



VofcTT 

2fc 



1 - 2a 
Vafc372 



+ ln(Vafc + VTTofc) 



fc3/2 



■ arctan 



(1 - a)k 
1 + afc 



(C3) 



/2"(a,fc) = 



Vl + afc(afc - 1 + 2a) 

2fc(fc + 1) 
V^(3 - 4a) 
2fc375 



ln(Vafc + vT+afc) 



7tSi 



VT^(1 - 4a) 
2fc372 



arctan 



(l-a)fc 
1 + afc 



(C4) 



M 



-S2 



(B12) 



By defining a new coordinate system in which the 
2;-axis is parallel to G, one find that the integrals of 
Eq. ljB12p involve gaussian integrals of the form 

y"dsexp(-s2) (\G\s,fQ{±s,)G,s, = ^IGpG, (B13) 

and 

j dsexp(-s2) {\G\s,fQ{±s,) = ^\G^\ (B14) 

which finally leads to Eq. H23() . The equation for rota- 
tional energy is derived following exactly the same pro- 
cedure. 



joi (a, k) =aif ( V^) + (1 - a)n(-fc, V^) (C5) 

where K{x) and T\{x^y) denote the complete elliptic in- 
tegral of the first kind and the incomplete elliptic integral 
of the third kind respectively. 



Jr{a,k) = 



K{ 



ak 
{l-\-ak) 



(i-«) 

2(fc + l); ^{l + ak) 



(l-a)v/(l + afc) 



ak 



2{k+l) 1 
(1 -a)(2afc + a + fc + 2) 
2(fc + l)V(l + a^) 



ak 



n(- 



1 + afc' 
(C6) 



10 



where E{x) denotes the complete ehiptic integral of the 
second kind. 



vT+ak / ak 
k W l + ak' 

I — a I k 

=n 



fc(fc + 1)^(1 + afc) yk+l'\ll + akj 

(C7) 



and finally 



J^'\a,k) = K{\ -) 

^ ' ' 2fc(fc + l)VTT^ Wl + afc^ 



\/l + ak{2ak + 3a — 1) , / afc , 
2fc(fc + 1) W 1 + afc^ 



2fc(fc+l)V(l + aA:) yk + V \ 1 + ak 

(C8) 
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